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1 Introduction 

The Korteweg-deVries equation is a simple mathematical model for gravity 
waves in water, but it fails to model such fundamental physical phenomena as 



the extreme wave of Stokes p3| . The failure of weakly nonlinear dispersive 
equations, such as the Korteweg-deVries equation, to model the observed 
breakdown of regularity in nature, is a prime motivation in the search for 



alternative models for nonlinear dispersive waves | 20| , [gi] 



In 1976 Green and Naghdi |13| derived a system of water wave equations 
to model fluid flows in thin domains, such as internal waves in coastal re- 
gions. The Green-Naghdi equations have a Hamiltonian structure, and in 
1993 Camassa and Holm |7J used scaling and an asymptotic expansion to 
obtain, in the one dimensional case, an approximate Hamiltonian which is 
formally integrable by the method of inverse scattering. The strongly non- 
linear equation they obtained, 

u t - -u xxt + ~(u 2 ) x - -(u 2 x ) x - -{uu xx ) x = 0, (1) 

supports solutions, dubbed "peakons", that are continuous but only piece- 
wise analytic. Equation (JJ) had originally been obtained by B. Fuchssteiner 
| 12| by the method of recursion operators in 1981. He showed the equation 



was Hamiltonian, but gave no physical interpretation, nor an isospectral op- 
erator; and the equation attracted no special attention until its rediscovery 
by Camassa and Holm. 

Motivated by the form of traveling wave solutions of (|1|), Camassa and 
Holm proposed solutions of the form (in the normalization used in this paper) 

1 n 
u{x,t) = ^^2 m j(t) exp(-2|x - Xj(t)\). (2) 

i=i 

to represent n interacting traveling waves. They substituted this Ansatz into 
(]l|) and obtained a Hamiltonian system of equations for rrij, Xj] the Hamilto- 
nian is obtained directly by substituting (fj) into the formal Hamiltonian for 
([]]). This system describes geodesic flow on a manifold with metric tensor 
g^ = exp(— 2\xj — Xj\). A solution (Q) can contain both peaks rrij > and 
antipeaks rrij < 0; at large positive or negative time it is asymptotic to a 
superposition of single travelling waves: peakons and antipeakons. 

In P the Hamiltonian system for two peakons is integrated, and explicit 
expressions for the relative position x\ — x<± and relative momentum pi — p 2 



are obtained. A qualitative analysis of the interaction of peakons, and of the 
collision of a perfectly antisymmetric peakon-antipeakon pair, together with 
a number of numerical studies are carried out, and a formula for the phase 
shifts of the interaction of two solitons was obtained. 
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Figure 1: A peakon moving from left to right collides at time to — -2310 
with an antipeakon moving right to left, computed from the explicit formulas 
using Matlab. Here A_ = — 1, A + = 2, a±(0) = .5. The slope becomes infinite 
at the instant of collision. Sharp results on the steepening of the slope are 
given in §[7|. The vertical scale is exaggerated; the peakons are n ot as sharp 
as they appear here. 



The strongly nonlinear dispersive wave equation ([]]) possesses a remark- 
able mathematical structure and is the subject of a steadily growing litera- 



ture. In || we used the theory of continued fractions and formulas of Stieltjes 
f22fl (based on prior results of Frobenius |TT|) to give algebraic formulas for 
the multipeakon solutions (all rrij > 0). These formulas led to explicit ex- 
pressions for the asymptotic behavior of the positions {xj} and momenta 
{rrij} as t — > ±00, and thus to explicit calculation of the scattering shifts 
undergone by the peaks as a result of their interaction. 

In this paper we give a brief but self-contained discussion of the spec- 
tral problem associated with ([[]), detailed proofs of the results announced 
in ||, and an extension of those results to the case of peakon/antipeakon 
interactions. In the remainder of this section we outline the paper, place it 
in the context of related work, and describe the conclusions concerning the 
solutions (0). 

In §0 we obtain the multipeakon/antipeakon solutions as a restriction of 
([]]) to a singular subclass. This necessitates a reinterpretation of the Lax pair 
in the sense of distributions. In §|3| we describe the Liouville transformation 
that transforms the spectral problem on the line to a formal density problem 
for a finite string, first investigated by M. G. Krein. In §^ we describe the 
Weyl function, which encodes the spectral data, and its expansion in contin- 
ued fractions. In §^ we describe the formulas by which the continued fraction 
expansion is recovered from the Laurent series for the Weyl function. In §^ 
we use the explicit formulas to obtain the asymptotic behavior of the mul- 
tipeakon solutions and the scattering shifts of the peakons and antipeakons 
that result from their interactions. 

A closed form of the multipeakon/antipeakon solution is given (Theorem 
|7.3|) in terms of the orthogonal polynomials associated with the moment 
problem. 

Stieltjes' original theory was restricted to positive weights, but the for- 
mulas are algebraic and extend to weights of mixed sign, so long as certain 
determinants do not vanish. In the present context, these determinants are 
functions of time. Collisions of peakons and antipeakons occur precisely 
when one or more of these determinants vanish. At such a point some of the 
weights rrij become infinite, but the solution itself remains bounded through- 
out the collision. Moreover, there can be no triple collisions: at a given time, 
collisions occur only in distinct peakon-antipeakon pairs: a peakon moving to 
the right encounters and passes an antipeakon moving to the left. We prove 
these facts in §|7|, using the explicit solutions together with the oscillation 
properties of the orthogonal polynomials of the associated moment problem. 



A primary reason for the interest in (|I|) is that it is an integrable model for 
the breakdown of regularity, a phenomenon not modeled by the Korteweg- 
deVries equation. A number of qualitative results of a general nature on the 
steepening of the slope at the instant of breakdown have been obtained by 



analytic methods |§, fllCj], j9fl; McKean, lfL7|1 , has given a sharp result for the 
general case. The exact results obtained from the Stieltjes formulas yield 
the same result (Theorem [7.2|) concerning the steepening of the slope as the 
instant of collision of a peakon and antipeakon is approached. 

In ||18| , fl9|l , J. Moser applied the theory of continued fractions to fi- 
nite Jacobi (a.k.a. Toda) flows. He obtained an explicit solution for 2x2 
and 3x3 Jacobi matrices, though he did not use the full results of Stielt- 
jes. He nevertheless deduced qualitative properties of the solution to the 
isospectral deformations of finite Jacobi matrices, and calculated the scatter- 
ing shifts. The relationship of Moser's work, as well as the classical formulas 
of Frobenius and Stieltjes, to the Riemann theta functions is discussed in an 



unpublished note of H. P. McKean |jl6 



Alber et al. 0,0 obtained a number of special solutions of the Camassa- 
Holm equation, including peakons, n— solitons, and quasiperiodic solutions, 
by inverting the transformation to action-angle variables. They attacked the 
Hamiltonian systems directly, without the use of inverse scattering theory. 
They assert that the general inversion problem can be solved in terms of Rie- 
mann theta functions; but explicit results, such as the phase shift undergone 
by interacting solitons, are obtained only in the case n = 2. 

Moser points out that the moment problem was the forerunner of modern 
spectral theory. Boris Levitan once remarked to one of the present authors 
that it was Gelfand's interest in a continuum analogue of the classical moment 
problem that led to their celebrated solution of the inverse spectral problem. 
B. Simon in [[H|] appears to close the circle, and returns to the original view 
of inverse spectral theory as a continuous analogue of the moment problem. 

The theorem of Stieltjes and the method of continued fractions was first 
applied to an inverse problem for ordinary differential equations by M. G. 
Krein |TJj|, (|T5J. The method was extended to a special class of fourth or- 
der equations by V. Barcilon |3J, which suggests that there are interesting 
isospectral deformations of such fourth order equations. 

The multipeakon/antipeakon solutions are intimately related to Jacobi 
matrices, but differ from the Toda flows in a fundamental way. The Toda 
flow describes the dynamics of a lattice of constant weights joined by restoring 



forces; the conjugate variables are the relative positions of the masses, and 
their momenta. In the multipeakon flow, the weights rrij vary in time and 
may take on both positive and negative values; the conjugate variables are the 
positions of the weights, and the weights themselves. We plan to investigate 
the relationship of the Jacobi flows to the moment problem in a future paper. 

2 The Camassa-Holm equation 

Equation ([I]) is obtained as the compatibility condition for an overdetermined 
system J7|; in the present normalization we take 

Lo(z)f = 0, ^ = A (z)f (3) 

L (z) = D 2 - z m(x, t) - 1, A (z) = I u(x, t)\ D + -u x (x, t) 

where i6R, z G C. 

Differentiating the first equation in (|3|) with respect to t and the sec- 
ond twice with respect to x and setting f xxt = ft xx , we obtain, after some 
calculations, 

m t + (um) x + mu x = 0, 2m = 4u — u xx . (4) 

We assume here that m and its derivatives vanish at ±oo. According to 
the second equation in (|j) the same is then true of u and u x . 

Conversely, suppose that the function m evolves according to ([|), with 
m and its derivatives vanishing at x — ±oo. This is also the compatibility 
condition for two modifications of ([3]): 

L o (z)p o = 0, ^=( i 4 ( z )-I) Vo (5) 

L (^ = 0, ?^=(A (z) + ^)i> . (6) 

These evolution equations are consistent with the asymptotic conditions 

<fo{x,z) ~ e x , x -> -oo; (7) 

ip (x,z) ~ e ~ x , x — >• +oo, (8) 



respectively. It follows that there are unique wave functions <p (x,t,z), 
i/j (x,t,z) that satisfy (|5|), (|7|), and @, (§), respectively. 
The remaining x-asymptotics of tpo are necessarily 

<fo(x,t, z) ~ b(t, z)e x + c(t, z) e~ x , x — >• +00. (9) 

It follows from the time evolution of ipo that 

db dc 2c ,_ . 

* = °- * = —■ (10) 

The eigenfunctions for this spectral problem are, by definition, the functions 
tp for those values z = \ v for which b(X v ) = 0. Therefore fliTf ) implies that 
the spectrum {A^} is invariant under the flow. The coupling constants are 
the values c u (t) = c(t, X u ); they are characterized by the relation 

y^oO, A„) = c u ?p (x, \„). (11) 

The evolution is given by specializing (]T0|) : 

2c 



^i/ 



(12) 

Thus far we have tacitly assumed that m is a continuous density. The 
interesting case for present purposes is when m is taken to be a discrete 
measure with weights rrij at locations xf 



m 



^ Wj^, ^1 < x 2 < ■ ■ ■ < x n . (13) 



The two equations in (|j) are readily interpreted in the sense of distributions. 
As we note below, the function u of (f|) arises in the same way. It can 



be calculated explicitly from (|T3"D and the second equation in (|j). In the 
normalization of the spectral operator L we have chosen, we obtain precisely 
the expression (§). We note that u x has jump discontinuities on the support 
of the singular measure m, so the meaning of mu x in the first equation (f|) is 
a priori ambiguous. We shall show, however, that the derivative D can be 
extended to piecewise smooth functions in such a way that all the operations 
extend to the discrete case without change. The operations which led from 
the overdetermined pair of equations (H) to the nonlinear system (H) therefore 



continue to hold in the discrete case. In particular, the evolution of coupling 
coefficients ([E|) applies to the multipeakon solutions as well. 
Suppose that a function / has the form 

f(x) = fj(x), Xj < x < x j+1 , (14) 

where each fj belongs to C°°(R) and we have set xo = — oo, x n+ \ = oo. 
We normalize at the jumps by taking the average of the limits from left and 
right: 

fix,) = /fe+) ; /fe -' = </<»,». 

If g is a second such function normalized in the same way then the distribution 
derivative D satisfies the Leibniz rule: D(fg) = fDg + gDf. Moreover, 
D(fm) = fDm. 

In the presence of smooth t dependence, 

/OM) = fjfat), Xj{t) < x < x j+ i(t), (15) 

/ as distribution has ^-derivative 

f = ^ ~ X>(^'+) - /(^-)1 *Ay- ( 16 ) 

Moreover 

3=1 

In addition, we have 

, n n 

D(um) —— y^ u(xj)rrijS Xi = '^2u(x j )m j D5 Xj ; 

3=1 3=1 

n 

u x m =y^ j (u x {x j ))m j 8 Xj . 

3=1 

The first equation of @j is therefore 

n 

^ rhjS x . - m j x j D5 Xj + {u x (x j ))m j 5 Xj + u(x j )m j D5 Xj = 0. 

3=1 

8 



Setting the coefficients of the independent distributions 5 X and D8 X . equal 
to zero, we obtain the Hamiltonian system 



where 



OH 
3 drrij 


= u(xj), rrij = 


dH 

dxj 


-{u x {xj))mj 


I(x,m) = 


/ \ u2 + 7 M ^) dx ~~ 


1 n 


m j m k e~ 2 ^ : >~ x ^ 



(17) 



1 n 1 r°° 

= ~J2 u ( x i) m j = ^ u(x)dm(x). (18) 

/ j= i A J-oo 

3 A Liouville transformation 

The flow of the amplitudes rrij and the locations Xj of the peaks can be 
computed explicitly, by virtue of the inverse scattering approach. As we shall 
see, there are n eigenvalues \j and coupling constants Cj. The eigenvalues 
are fixed under the evolution and the coupling constants evolve according to 
(|T2|). In order to use this information to recover the time dependence of rrij 
and Xj we find it convenient first to transform the spectral problem. 

In we obtained a Liouville transformation which transforms the spec- 
tral problem L (z)ip = to a "density" problem for a finite string, in which 
the density can take on both positive and negative values. The density prob- 
lem with positive sign has been studied extensively by Krein |14| , ]15| . In the 



research announcement || we used the classical results on continued fractions 
to obtain solutions of (0) in the case m is a discrete measure, with weights 
of one sign. The argument carries over to weights of arbitrary sign. 

For completeness, we repeat the details here. The first step is to transform 
(a multiple of) the operator D 2 — 1 that occurs in (^) into the operator D 2 on 
the interval (—1, 1); then the wave functions for discrete m become piecewise 
linear. The appropriate coordinate transformation is 

y = tanh x, — = p(x) = — - 1 — = 1 - y 2 . (19) 

dx cosh x 

The operator 

y 2 L (z) = cosh 4 x {^-zm-l) (20) 



is selfadjoint in L 2 (M, p 2 dx). The operator 

U :L 2 ((-l } l) } dy)^L 2 (R } p 2 dx), (21) 

[Uf\{x) = p(x)- l ' 2 f(timh x) = (1 - y 2 ) 1/2 f(y) 
is unitary and carries (120) to 



L{z) = -^ - z(j(y), g(y) = m^anh" 1 y)(l - y 2 ) 2 = ^, (22) 

with Dirichlet boundary conditions. As noted above, the resulting spectral 
problem 

^(y)=zg(y)v(y), -1< y < 1; v(-l) = = v(l), (23) 

is familiar when g is positive: it determines the natural vibration frequencies 
a string with mass density g. 

The Liouville transformation takes the wave functions (fo, ipo to the cor- 
responding functions for the problem 



V /'-zg<p = 0; p(-M) = 0, ^'(-1,0) = 1; (24) 

ij/'-zgil> = 0\ </>(M)=0, ij>'(l,z) = -l; 

here the primes denote derivatives with respect to y G [—1,1]. The eigen- 
values {Aj,} and coupling constants {c u } are the same as for the original 
problem. 

In particular, the evolution of the coupling coefficients is preserved under 
the Liouville transformation. 

The discrete measure (|13D is transformed into the discrete measure g on 
the interval (—1, 1) given by 

n 
9 = ^2gjd yj , "I = Vo < Vl < ■ • • < Vn < Vn+l = 1- (25) 

The terms in ([13]) are related to the terms here by || 

m j = g j (l-y 2 ), *, = ± tog (I±l£) . (26) 

10 



We remark that the Liouville transformation is a canonical transforma- 
tion from the equations on the line to a Hamiltonian system on the interval 
(—1,1), but we shall not need that here. The multipeakon solution is ob- 
tained by solving the inverse spectral problem for the discrete string, and 
then pulling the solution back to the real line, using fl26f) . 

In the next two sections we treat the direct and inverse spectral problems 
for the discrete string. 

4 The Weyl function 

Equation (p3|) implies that <p is continous and piecewise linear, with jump dis- 
continuities in the derivative at the points {%}. Then both sides of ( p3[ ) are 
interpreted in the sense of distributions. We denote left and right derivatives 
with respect to y by D±: 



D+viVj 



<p(yj+i) ~ <p(Vj) 



h 



D-<p(yj 



<p(yj) - y(yj-i) 



where I 



Vj+i ~ Vji < j < n. Then (p3|) becomes 

D+cpiVj) ~ D-<p(Vj) = z 9j ¥>(%•), l<j<n. (27) 

We note that the spectral problem may be written in the matrix form 



Jq = zGq, 
where q = (q u . . .,q n )*, qj(t) = tp(vj(t),t), 
/6i oi ... \ 



(28) 



J 



Oi b 2 a 2 
a 2 h 







dn-1 b n J 



G 



fgx 







0\ 





92 








V° 







9n,J 



and Oj = 1/lj, bj = — (aj_i + %). The spectral problem for the multipeakon 
solutions thus differs from that for the Jacobi flows, which simply takes the 
form Jq = zq. 

The wave functions tp(-,z) may be constructed recursively. We fix z and 
set 

q j = ^(y j ,z), Pj = D-(p(y j ) } 



11 



so that 

1j ~ Qj-i = Pjlj-u Pj ~ Pj-i = z 9j-iQj-i- ( 29 ) 

We begin with q = and pi = 1. Then qj and pj are polynomials of degree 
j — 1 in 2. In particular, (p(l, z) = q n+ i and D_y?(l, 0) = p n+ i are polynomials 
in 2 of degree n. 

Theorem 4.1 T7ie set of roots of<p(l, z) is the spectrum of the problem (p3|) 7 
equivalently (p8|) . T7ie roote are rea/, simple, and non-zero. The number 
of positive (resp. negative) roots of ip(l,z) equals the number of negative 
(resp. positive) terms among the gj. 

Proof: The eigenvalues are given by the zeroes of (p(l,z). Any other non- 
zero solution of (|23|) which vanishes at —1 is a scalar multiple of ip(y,z). 
Therefore there is at most one eigenfunction for each value of z, and the 
geometric multiplicity is one. Moreover, cannot be an eigenvalue, since 
<p(y, 0) = 1 + y. The eigenvalues are real since J and G are real, symmetric 
matrices and J is negative definite. 



Consider (28) with G — I. This is a pure eigenvalue problem for the 
symmetric tridiagonal matrix J, so it has real spectrum. It corresponds to 
our spectral problem with all gj — 1; so if z > 0, then the slope of (p(-, z) is 
non-decreasing from left to right and necessarily ip(l,z) > 0. Consequently 
J is negative definite and has the form — B 2 , where B is positive definite. 
Then Q28D is equivalent to B~ l G B~ 1 w = —z^ 1 w, w = Bq. Since B^ 1 GB^ 1 
is symmetric, the spectrum is real and the algebraic multiplicity of each 
eigenvalue is its geometric multiplicity, which we have shown to be 1. 

To obtain the conclusion about the signs of the roots, we set B(s) = 
(1 — s)I + sB, < s < 1. Each B(s) is positive definite, so by the previous 
argument the eigenvalues of B(s)~ 1 GB(s)~ 1 are simple and non-zero. There- 
fore the number of positive roots is independent of s. Comparing s = 1 and 
s = 0, we obtain the result. q.e.d. 

We denote the roots of <p(l, z) by {\j}- Note that <p(y, 0) = 1 + y, so 
</?(l, 0) = 2 and 

n 

¥>(M)=2n(l-£). (30) 



12 



The Weyl function associated with (|27D is 



w(z) 


-. n n 


z 


2z 4-i z- \j ^ z - A 7 




ao = ~, A = 0. 



D_^(l,z) p n+ i 

w z = — n — r - = — • 31 

<p(l,z) q n +i 

It will be more convenient to work with the modified Weyl function w(z)/z; 
in particular we use its partial fractions decomposition 



(32) 

Z LZ ' ' Z — Aj ' ' Z — Aj 

where we have set 

a,n = 

2' 

The scattering data for the spectral problem (p7|) consists of the spectrum 
{Aj} and the coupling constants {cj}. Recall that the coupling constants 
relate the eigenfunctions ip to their counterparts if) normalized at y = 1 by 
V>(1, z) = 0, D_V>(1, z) = -1. In fact 

<p{y,><j) = Cjil>{y,\j). (33) 

The Wronskian of two solutions of Q2"7| ) is constant in the intervals y 3 - < 
2/ < 2/j+i an d continuous across the yj. Evaluating the Wronskian W(ip, ip) at 
y = ±1 establishes that f(l,z) = ip(—l,z). Differentiating (|33|) with respect 
to y and setting y = 1, we obtain 

c i = V(1,A,). (34) 

The residues %, j > 1, in ( |3"2"D are determined by the scattering data. 
Combining (|34]) and (0), we obtain 



Oj 



r^rrrr = IIId-VA*)- 1 . (35) 

Aj ^(l,Aj) 2 " 



Thus under the nonlinear evolution given by ([[]), 

2a 



0;' = -^, aj(t) = a i (0)e" 2 * /Aj . (36) 

Ai 
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Theorem 4.2 The residues clj in the partial fractions decomposition ( |3"2"D of 
w(z)/z are positive and satisfy 



X 2 V 2 z (l,X J )a j = J y 2 (y } X 3 )dy. (37) 

Proof: Differentiation with respect to z commutes with D, the distribution 
derivative, so 

D 2 (p = zipg ) D 2 (p z = ZLf z g + Lpg. 

We multiply the first of these equations by <p z and the second by <p, subtract, 
and integrate with respect to y. When z = Xj we obtain 

ip'(l, Xj) (p z (l, Xj) = - ip 2 (y, Xj) dg(y), (38) 



where dg is the measure defined by the distribution g defined by (25). Inte- 
grating the equation ipD 2 {p = Xjip 2 g with respect to y, we obtain 



i pi 



V9 /2 (y, A,) dy = -A, J <p 2 (y, Xj)dg(y). (39) 

Combining (J35|) , Q38|) , and (|3~9"|), we obtain (|37|). q.e.d. 



5 A theorem of Stieltjes 

For the inverse spectral problem, we assume that the eigenvalues {Xj} and 
the residues {a,} are known, and seek to determine the constants {gj} and 
the points {%•}; in place of the latter we may look for the subinterval lengths 
lj = yj + i — yj. A first step in the process is to determine the Laurent 
expansion at infinity of the modified Weyl function. This is easily obtained 
from (|32"D by expanding 



"3 



z — A,- *-^ z 



\i 



yk+l ' 



fc = 



The result is the following. 

Lemma 5.1 The modified Weyl function has the Laurent expansion 

w(z) ^ (-l)M t 



E^P- («» 



z 

k=0 
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where 



A k = E(-^) fe %- ( 41 ) 

3=0 



The Weyl function itself can also be written as a continued fraction, pA 

Lemma 5.2 The Weyl function is 

w(z) = l — . (42) 

k + : 



zg n 



ln-1 + ' 


1 

' + i 

L 





Proof: We use (^9]) inductively: 

Pi = 1, q 1 = l , q 2 = q± + hp 2 , P2 = Pi + zg x l Q ] 



P2 P2 1 



Q2 k + P2h , 1 

h+ 



1 
zgi + y 



Assuming (H) for {vj}j <n and {gj}j <n , we adjoin y n > y n _i and g n . Then 
Pn+i =Pn + zg n q n , q n+ i = q n + l n p n +i] and 



Pn+l Pn+1 



1 



<?n+l Qn + InPn+l q n q-n 

l n "T in ~T 



Pn+1 Pn + ^n^n 



1 



(43) 



Pn 
?n 
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Hence fl4"2"|) follows by induction, q. e. d. 

Dividing by z gives the continued fraction decomposition of the modified 
Weyl function: 



w(z) 



(44) 



Lz + 



9n + 



1 



J„_i 2 H h 



1 



01 



1 
Z £ 



A classical result of Stieltjes |22| recovers the coefficients of the continued 
fraction ([E|) from the Laurent expansion oiw(z)/z at infinity 

A key role is played by the Hankel matrix whose entries are the quantities 
Ak defined by (|4l|). In the original theory of Stieltjes, the Ak are all positive. 
This corresponds to the pure peakon case, in which all the weights gj are 
positive and the eigenvalues Xj are negative. In the peakon- ant ipeakon case 
the weights are of both signs, and the Ak are now of both signs; nevertheless, 
- and this will be an important point in the analysis -- they are still the 



moments of a positive measure. In fact (|4T| ) may be rewritten as 

nfty\ k d{i(\), (i = ^2a j 6-\ j , (45) 

i=o 

where again we take A = and a = 1/2. 

In view of fl4"T| ) and (^) , we see that the A,- are rational functions of the 
scattering data. The following result provides explicit formulas. 

Theorem 5.3 [Stieltjes] The Laurent series ( f4"UD can be uniquely developed 
in a continued fraction 

1 



where 



>2k 



1 
b lz + — 

b 3 z + . 


• 




(A°J 2 , 

AlAl ' °2fc+l " 


(AD 2 
A° fe A° fc+1 



(46) 
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Moreover 



61 + h + ■ ■ ■ + b 



2k+l 



A° fc+ i 



(47) 



Here the A\ are certain k x k minors of the infinite Hankel matrix 

(Aq A x A 2 A 3 ...\ 

A x A 2 A 3 Ai .. 

A 2 A 3 A 4 A 5 .. 

A 3 A 4 A 5 A 6 .. 



H 



\- 



J 



A l k is the determinant of the k x k submatrix of H whose (i,j) entry is 
Ai +i+ j^ 2 - By convention, A° = 1. 



Equations (HBI) and ( f£7| ) appear in [EJ|, equations (7) and (11), respec- 
tively. By comparing the continued fraction in ( [P]) with that in Theorem 
we obtain 



A° A 



0<j<n; 



(48) 



(A°- J+ i) 2 

9j Al Al ' 

71— j'+l n— j 

while from (||7|) and (H3) we obtain 

j/,- = 1 - (/„ + /„_i H \-lj 



l<j<n; 



A 2 ■ 
A 



(49) 



(50) 



These results allow us to characterize the range of the forward spectral 
map in terms of conditions on the A&. 

Lemma 5.4 For even I, the minors A l k are strictly positive, < k < n, as 
is A° +1 . For odd I, and < k < n, if all the weights gj have the same sign, 
then A k is non-zero and has the opposite sign, while A l n+1 = 0. 
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Proof: The quadratic form associated with the k x k submatrix with upper 
left hand element A\ is 



fe-i „,• fe-i 



i,j=0 J ~°° i,j=0 J ~°° 

where 

fe-i 

j=0 

Thus the form is positive definite when / is even and k < n: a non-zero 
polynomial of degree k — 1 cannot vanish at the n non-zero points — Xj in 
the support of dfi. (When I = the point A = must also be considered.) If 
all these points — Xj are positive (respectively negative) then the form is also 
positive (respectively negative) for odd I. The associated determinant is A l k . 
Finally, for odd I and k > n, there is a nonzero polynomial Yl c j^j °f 
degree k — 1 that vanishes at each Ai, ... ,A n , and it follows that the vector 
(c , . . . , Cfc_i)' is in the null space of the matrix associated with A l k . q.e.d. 

In the pure peakon or antipeakon case, i.e. when all the weights rrij have 
the same sign, Theorem [|J] says that all eigenvalues Aj have the same sign. 
It follows immediately from Lemma |5.4| that the minors A] cannot vanish; 



hence lj > for all time, and there are no exact collisions: the distances 
Xk+i — Xk are always strictly positive. 

Theorem 5.5 The real non-zero constants {Xj,Cj}, j — 1, ...n are scatter- 
ing data for (|23D if and only if the Xj are distinct, the aj given by (|35D are 
positive, and the moment matrix constructed from the associated measure has 
the property that the determinants A 1 -, < j < n, do not vanish. 

Proof: Suppose that the {A^}, {cj} are scattering data. We showed in §[|that 
the eigenvalues are distinct and the aj are positive. It follows from Lemma 
|x4| that A° > 0, < j < n+ 1. Since the lj are positive, it follows from (|48l) 
that A} ^ 0, < j < n. 

Conversely, starting from distinct {A^} and {cj} such that {aj} are pos- 
itive and Aj ^ 0, (^) and ( |4^ ) can be used to define {lj} and {<?•,}• The 
Laurent series ( |40D corresponds to the continued fraction (0). The proof of 
Lemma ( |5.2| ) is reversible The associated continued fraction P2| ) terminates, 
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because A* +1 = 0, and gives the Weyl function for the spectral problem with 
constants {gj} and {lj}. q.e.d. 

To relate the measure m to the scattering data, we begin with some no- 
tation and an identity. We denote by Ak and A\ the corresponding moments 
and determinants with respect to the modified measure 

n 

Ji = ^2a j 6-x i - (51) 

3=1 

These coincide with the A k and A^ except that A = 1/2 + A , and conse- 
quently 

A2 = A2 + iAti, k>l. (52) 

It follows from (|0|) and (0) that 

A 2 2A° 

1- Vj = -^-, l+^^^P 11 - (53) 



Combining (^) with (|U) and ( |2T)| ) we obtain the following. 

Theorem 5.6 The weights rrij and positions Xj associated with the distribu- 
tion m in fll3|) are given by 

i /9A° s 9A° A 2 



o "° \ A2 /' -'-J Al Al 

We illustrate with the cases n = 1, 2, 3. For n = lwe have 

1 /2A°\ 1 



xi= 2 log iArJ = 2 log2ai 

2A? Aj 2ai 2 



mi " AJ AJ " -A l0l " A! 
For n = 2 the result is 

1 2(Ai - A 2 ) 2 aia 2 1 

2 g A 2 a! + A 2 a 2 ' X2 ~ 2 



x 1 =-log— ^2 —2 , x 2 = -log2(ai + a 2 ) 



2(A?Oi + A 2 a 2 ) 2(oi + a 2 

m i = - t— r~n. — ; r> m 2 



AiA 2 (AiOi + A 2 a 2 )' ~ Aiai + A 2 a 2 
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Finally, for n = 3 we have 



X-2 



x __l lo 2(Ai - A 2 ) 2 (A 2 - A 3 ) 2 (A 3 - A 1 ) 2 a 1 a 2 a 3 
= 2 lQ g( ^" )> x 3 = -log(2^a i ); 



mi 



A1A2A3 Ylj<k ^j^ki^j ~ Afc) a^afc 



2 E A'% Ej-< fc (Ai - Afe) 2 ^-^ 2 £ a, 

m 2 = ~^TT F^ , , 7T r^ > m 3 - 



Moser Jl9| applied the theory of continued fractions for Jacobi matrices 
and obtained explicit solutions similar to these for the isospectral flow of 
3x3 Jacobi matrices in a special case. 

6 Asymptotics 

The function u in (H) is a superposition of single peakons (m^ > 0) and an- 
tipeakons (m& < 0) with peaks and troughs at the Xk] however both the rrik 
and the Xk vary in time. We show now that u is asymptotic at large times 
to superpositions of non-interacting peakons and antipeakons with constant 
heights/depths: At large negative time peakons are found to the far left, 
decreasing in height from left to right, travelling to the right at speeds pro- 
portional to their heights, with antipeakons to the far right, increasing in 
depth from left to right, travelling to the left at speeds proportional to their 
depths. At large positive times antipeakons are to the left, with decreasing 
depths, and peakons to the right, with increasing heights. The difference, 
relative to a pure superposition of such travelling solutions, is that each peak 
or trough has undergone a phase shift. The total (asymptotic) phase shift 
for each peak or trough is precisely the sum of the phase shifts that would 
occur in interactions with each of the other terms separately (in a calculation 
with n — 2); see fl62|) . 

The key to the analysis of the long-term asymptotics is the evaluation of 
the determinants A l k and A°. 
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Theorem 6.1 The determinants A l k of the moment matrix for the measure 
d\i in ( f45|) are given by 

A[= J2 « J (-A)^j, (55) 

Jc{0,l,...,n},\J\=k 



where 

« J = II °J' AJ = II A J' nj = 1 1 ( X i ~ Xm 

jeJ je.J j,m€J,j<m 



The determinants A° are 

ai= Yl aJ ^- ( 56 ) 

JC{l,...,n}, \.J\=k 

The proof of this theorem is given at the end of this section. We begin 
by using it to deduce the asymptotics, We number the eigenvalues Xj so that 

A m+ i < • • ■ < A n < A = < Ai < A 2 < ■ ■ ■ < A m . (57) 



An examination of the formulas Q55D and (52) in the light of (^) and the 



assumption (|57|) yields the following asymptotics of the determinants. 
Lemma 6.2 As t — > — oo, if I > and k < n, then 



k k 

Ai~(n«i(°)(-*i)')( II (Ai-A r ) 2 )ex P (-^ A7 

3=1 l<j<r<k j=l 3 



2t\ 



If k < m then 



k 



2t\ 

7t A J 



A2~(n«i(0))( II (A,-A r ) 2 )ex P (-^£). 

3=1 l<j<r<k j= 

Ifm<k<n+1, then 

fe-i fc-i 



A °-^(II%(0))( II (A,-A r ) 2 )exp(-^ 



2t 



A,- 

j=l l<j<r<fc j=l J 

In a// cases, for k <n 

k k 



A°~(n%(0))( II (A, - A r ) 2 ) exp ( - ^ ^' 

1=1 l<7<r<fc 1=1 ■?' 
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Proof: We discuss only the first case; the remaining are treated similarly. The 
sum in (|55|) is taken over all subsets of integers J of size k. The exponential 
term that corresponds to any such subset is 



""("££) 



Since —t — > oo, the dominant term is that for which J = {1, . . . , k}. Hence 

k 



A&) -n^-wH') 1 n (a s -aj 2 

j=l s<m<k 

k / k 

= na i (0)(-A i )' J] (A s -A m ) 2 exp -£^ 

j=l s<m<k \ j=l ■* 

and the first result follows. 

Lemma 6.3 ^4s £ — > +oo ; if I > then for k < n 

Ai~( II %(0)(-^)( II (A, - A r ) 2 ) exp ( - £ |). 

j>n—k n—k<j<r j>n—k J 



If k < n — m then 



n 



K-(U o,(0)) ( II (A, - A r ) 2 ) exp ( - £ £ 

j>n—k n—k<j<r j>n—k ^ 

If n — m < k < n + 1, then 

*i~\( n ww n (a, - a,) 2 ) ex P ( - £ f> 

j>n— fc+1 n-t+l<i<f j>n— fc+1 ^ 

In a// cases, for k < n 

n n n _ 

A2~( II «i(°»( II (A, - A r ) 2 ) exp ( - £ £' 

j>n—k n—k<j<r j>n—k ^ 
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Combining these results with (|5.6|), we obtain the asymptotics of the problem 
on the line. 



Theorem 6.4 As t 



-oo, 



x n _ i+ i(£) ~- — + -lo 



A, 2 



20,(0) n (1-^ 

fc=l 



k 



m .(i) ~- — 



(58) 
(59) 



As t — > +oo ; 



Xj(t) 



H + 2* 



n \ 2 

2%(o)n(i-£ 

fe=j+i 



m n _ i+ i(t) ~ - — . 
A i 



(60) 
(61) 



It follows that the solution u(x, t) is asymptotically a sum of free peakons 
(Xj < 0) and antipeakons (Aj > 0). The term with asymptotic velocity 
— 2/Aj undergoes a phase shift due to the interactions in the amount 



fe=i+i 



^1 



i-i 



£>g 1 



fc=i 



A/,. 



(62) 



Proof of Theorem 153]. The proof is based on the following lemma. 



Lemma 6.5 Let A^, k = 0,1,2, ... , be the moments of the discrete measure 



A^J^A 



i=o 



where the weights bj are non-zero, but need not be positive. Then the deter- 
minant 



A 



A A x 
A 1 A 2 



A), 



A n +i 

A n A n+ i . . . A2 n 



(63) 
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is given by 

n n 

i=0 k>j 

Proof: A° +1 is a polynomial in the bj and Xj, so it is enough to verify the 
identity (6.10) when all the bj are positive. Under this asssumption, A° +1 is 
a moment matrix for a positive measure with weights at the Xj. By Lemma 
|5.4| , A° +1 is positive if and only the Xj are distinct. Therefore each root 
Xj — Xk is double, and A° +1 is divisible by EL^G^ — ^k) 2 - Comparing the 
total degree in the Xj, we find 

A° +1 = c n+1 (6 , • • ■ , K) ]J(X k - A,) 2 . (65) 

k>j 

The coefficient is a polynomial of total degree n + 1 in the bj and vanishes 
if any bj = 0, hence has the form c n+ in&.,-. The constant c n+ i may be 
determined inductively: we take A n = 0, so only the moment A has a b n 
term and A° +1 is the product of b n and the n x n minor A 2 . The latter is 
A° for the measure 

n-l 
3=0 

and it follows that c n+ i = c n — ■ ■ ■ — c\ — 1. q.e.d. 

Note that any larger moment matrix for the measure d[i supported at n+1 
sites vanishes, e.g. by introducing additional sites with weights b n+ k = 0. 

The formula (|55|) for / = 0, k = n + 1 coincides with (|6~4|), under the 
replacement Aj i— » — Aj. The remaining formulas may be derived from this 
one. First, suppose that / = and k < n. We express the determinant A° as 
a sum of terms indexed by the J's; each term is itself the determinant of the 
k x k matrix associated with the measure (53). This proves ([55]) when / = 0. 



When / > 0, A^, coincides with a determinant of the form A° taken with 
respect to the measure 

n 

^2(-Xj) l aj8. Xj . 

j=0 

Therefore the remaining formulas of (|55|) are consequences of the formulas 
for / = 0. q.e.d. 
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7 Collisions and close encounters 



The phase shift formula fl62|) suggests that interactions occur in pairs. This 
is borne out by an analysis of close encounters (xk+i(to) — Xk(ta) positive but 
locally minimal) and collisions (xk+i{to) = Xk(to))- When only peakons or 
only antipeakons are present, collisions cannot occur: the determinants A] 
do not vanish; see Lemma 5A. The equation ij = u(xj) in (|I7]) shows that 
an overtaking peakon must be higher, and that after the event, the higher 
peakon must be to the right; whereas an overtaking trough must be lower, 
and moving to the left. 




Figure 2: A close encounter, computed from the exact formulas; Ai = —1, 
A 2 = —2, a>i(0) = 1; a 2 (0) = .5., t = 0. After the encounter, for t > t , the 
peakon on the right picks up speed and moves away. 

When both peaks and troughs are present, collisions occur. At a collision, 
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some A*_ fc = and the solution of the inverse problem for (|2"3| ) breaks down: 
the terms m^ and mfc +1 blow up. We have Aq = 1, so m n cannot become 
singular unless m n _i does also, and (|55| ) implies that A^ > 0, so mi cannot 
become singular unless 1712 does also. 

We show by a direct analysis of the inverse problem that u has a well- 
behaved continuation throughout the collision. (This can also be seen from 
the fact that the Hamiltonian ([0|) is conserved under the flow, so that the 
W 1 ' 2 norm of u, and hence the L°° norm, remain bounded.) 

The qualitative facts about collisions between peakons and antipeakons 
are summarized in the following theorem. 

Theorem 7.1 Suppose that some of the determinants A^_ fc vanish at t = 
to. As t approaches to, collisions occur within distinct peakon-antipeakon 
pairs. For t < to each such pair consists of a peakon with peak at Xk and an 
antipeakon with trough at Xk+il for t < t the trough is at Xk and the peak is 
at Xk+l- 

As t — ► to± the function u(x,t) converges uniformly (with respect to x) 
to a function that has the same form (but with the Xj no longer distinct). 

These aspects of the peakon-antipeakon collision are illustrated in the 
figure in the introduction. 

Although u(x, t) is well-behaved in the supremum norm at collisions, this 
is not the case for the derivative. 



Theorem 7.2 Suppose that A n _fc(£o 

O: 



u x (x,t) 



t-to 



0. Then on the interval {xk,Xk+i), 
+ 0(1) as t^t , a>0. (66) 



The proofs of Theorems |7.1| and [Z.2| use some facts from the classical 
moment problem; proofs may be found in the monograph by Akhiezer, |[|, 
Chapter 1. Given the positive measure d\i of (|45|) , one obtains a sequence of 
n + 1 polynomials {Pj(X), < j < n}, such that Pj has degree j and the Pj 
are ortho normal in L 2 (M, dfi). Explicitly 



P,(A) 



(aJhAJ) 



1/2 



Ao 
At 



A, 

A 2 



Ai 
A j+ i 



A.i-i Aj 



A 



2j-l 

A' 



< j < n. 



(67) 
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In particular, 

hence, by ([48]), a collision between the k th and (k + l) st places occurs precisely 



P j (0) = (-iy —- L— 0<j<n. (68) 



when the constant term in the polynomial P n _fc(A) vanishes. 
The -Pj(A) satisfy a second order recursion relation 

XPj(X) = &;P/+i(A) + djPjiX) + b^Pj^iX), l<j<n-l, (69) 

where 

(^+ 2 ) 1/2 



bj= V J A r >0, 0<j<n-l. 



The recursion relation implies a well-known formula of Darboux-Christoffel: 



8=0 

We shall use the limiting form 

i 
6, (j^P, - P i+ iP0 =J2 P ^ < j < n - 1. (70) 

8=0 

where primes denote differentiation with respect to A. 

Theorem 7.3 27ie multipeakon solutions may be expressed in terms of the 
orthogonal polynomials Pj{0, t) as follows: 

h = P <" m2 > ^- ^Jo)p H( o) (71) 

ft = 1 + <>»-j (^-j(0)P„- J+ i(0) - P„- 3 (0)P^- +1 (0)). (72) 

Proof: The formulas follow immediately from (f48|) , (49), Q68[) and (ffO|). 



An immediate consequence of (|70|) and the assumption that -Po(O) ^ is 
that no two consecutive Pj{0) can vanish simultaneously Therefore no two 
consecutive A] can vanish, < j < n — 1; and collisions can occur only in 
distinct pairs rrik, rrik+i- Moreover, if Pj(0) = then Pj(0) 7^ 0. From now 
on we consider the time-dependent case, but we shall mainly suppress the 
time dependence in the notation, e.g. Pi = p(A,£), P(0) = Pj(0,£). 
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Lemma 7.4 Suppose that the measure ( f£5D evolves according to (|3T)|) ; anc? Ze£ 
{P,(A,£)} 6e the associated orthogonal polynomials. Suppose that Pj(0,to) = 
0. Then the derivative Pj at t = to satisfies 

P J + 2P J /\ = P;{0)Y f P i (p)P t . (73) 

i<j 

Proof: Let ( , ) denote the inner product with respect to dfi. Differentiating 
(Pj,Pi) = 5ij, with respect to t (and suppressing the time variable) gives 

= (Pj, Pi) + (Pj, Pi) + J2 2 ^7 Am) ^i(-Am) a m . (74) 

We assume that Pj(0) vanishes at t — t ; then Pj/X is a polynomial of degree 
less than j with constant term P'(0) and we can add a summand at A = 
in (|74|) to obtain 

(P„ P) + (P i; P) + 2(P i /A, P) = Pj(0)P(0). 

Note that (Pj, Pj) = if i < j. It follows that the polynomial Pj + 2Pj/A has 
degree < j and its expansion in the orthonormal basis {p} is (|73|). g.e.c?. 

Note that (|73|) implies 

p,(o) + 2p;(o) = p;(o)]Tp(o) 2 . (75) 

i<j 

The next lemma implies that the zeros of the Aj(t) are simple. 



Lemma 7.5 Under the assumptions of Lemma 7.4, Pj(0)P,'(0) < ait = to- 



Proof: Since the polynomial Pj + 2Pj/A has degree < n, it cannot vanish 
at all the non-zero points — A m in the support of the measure. Therefore we 
have a strict inequality 

i(P,(0) + 2P;(0)) 2 <||P J + 2A- 1 P J || 2 = ^P;(0) 2 p(0) 2 , (76) 

i<j 



where || || 2 is the norm on L 2 (R, dfi). The equality in ([T6|) follows from (|73|) 
and the orthonormality of the polynomials Pj. 
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It follows from ([7§) and (|7§) that 

-2P,(0)P,'(0) >P,(0) 2 . q.e.d. 

The next lemma allows us to determine the direction of the sign change 
for interacting peakon/antipeakon pairs. 



Lemma 7.6 Under the assumptions of Lemma VTA at t = to 



|(P J+1 (0,t)P,(0,t)) > 0, ^(^(0,0^-1(0,*)) < 0, (77) 

for 1 < j < n — 1 and 1 < j < n respectively. 
Proof: By Lemma [7.5| and the assumption that Pj(0, to) — 0, the derivatives 



in (|77D have the same sign as 

-p i+1 (o,t )p;(o,to), -p;(o,t )p,-i(o,to), 

respectively. The conclusion ([FT]) follows from ([70|) . g. e.d 

The final lemma gives an algebraic proof that the singularities in success- 
sive gkS, hence in successive m^'s, cancel each other in a collision. 



Lemma 7.7 Under the assumptions of Lemma \7.$ the sum gu + gu+i has a 
finite limit as t — > to- 



Proof: Let j = n — k. From (ffll) and (|69|) we obtain 

6 i P i+1 (0) + 6 i _ 1 P,_ 1 (0) dj 



6 i 6,_ 1 P i+1 (0)P ? (0)P ? _ 1 (0) 6 i 6,_ 1 P ;+1 (0)P ; _ 1 (0) 



fl'fc + S'fc+i 



The denominator does not vanish for t near to, and it follows from fl69|) that 
dj = (XPj,Pj), which has a finite limit, q.e.d. 



Proof of Theorem \7.1\ As noted above, ( |7TJ|) implies that collisions occur 
only in pairs m^, m,k+i, where A^_ fc vanishes. According to Theorem |5.6| 
and equation (p5|), m^ and m^+i have the same signs as — P /+ i(0)P ? (0) and 



— Pj(0)Pj_i(0) respectively, j = n — k. It follows from (|77|) that m^ changes 
from positive to negative and rrik+i changes from negative to positive. In 
other words, a peakon on the left changes places with an antipeakon on the 
right. 
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The assertion about uniform convergence of u(x,t) follows from Lemma 
[777| . In fact suppose that A*_ fe (£ ) = 0, so that m k and m fc+1 become singular. 
At the same time, however, y k — y k +i — * 0, so Lemma \T7\ and (|26| ) imply 
that the sum m k + m k+ \ has a finite limit. For such a pair 

m k {t)e- 2 \ x - Xk{t) \ + m fc+1 (£)e- 2|l -* fc + l(i)l 

= {m fe (t) +m fc+1 (t)}e- 2 l ;c -^WI +m fc+1 (t){e- 2|a; -^+ l(t)l -e-^-^W 1 }. 

The difference between the exponentials is uniformly 0(jx k (t) — x k+ i(t)\), so 
the sum converges uniformly to 



m fc (£ )e~ 2| ™ (to)l + m k+1 (t )e- 2lx - Xk +^ to)l , 



where 



m k (t ) =m k+1 (t ) = - lim{m fc (t) + m k+1 (t)}, 

L t—fto 

x k (t ) = \imx k (t) = lim x k+1 (t) = x k+1 (t ). 

t— >to t—*to 

The same argument applies to any other pair that become singular at t — to, 
so u(x, t) converges uniformly. q.e.d. 

Proof of Theorem |7.4 It follows from the results of this section that if 
K-M = 0, then 

x k+1 -x k = 0(4) = a {t - t ) 2 + 0((£ - to) 3 ); 

m k = --^- + 0(l); 
t — to 
0/.\ 

m k+1 = - — — + 0(1), 
t — to 

for some positive constants <x, which can be computed explicitly, from these 
facts and (0) we deduce that on the interval in question we may compute the 
derivative to order 0(1) by dropping all but two summands and replacing 
them by 

ft l /2x-2x k+1 _ ^x k ~2x\ 

t-to K ' 

whose derivative on the interval (x k ,x k+ i) is 4ai/(t — t ) + 0(1). q.e.d. 
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